Dynamical turbulent flow on the Galton board with friction and magnetic field 
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We study numerically and analytically the dynamics of charged particles on the Galton board, a 
regular lattice of disc scatters, in the presence of constant external force, magnetic field and friction. 
It is shown that under certain conditions friction leads to the appearance of a strange chaotic 
attractor. In this regime the average velocity and direction of particle flow can be effectively affected 
by electric and magnetic fields. We discuss the applications of these results to the charge transport 
in antidot superlattices and stream of suspended particles in a viscous flow through scatters. 
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It is well known that dissipation can lead to the appear- 
ance of strange chaotic attractors in nonlinear nonau- 
tonomous dynamical systems In this case the en- 

ergy dissipation is compensated by an external energy 
flow so that stationary chaotic oscillations set in on the 
attractor. Such an energy flow is absent in the Hamilto- 
nian conservative systems and therefore the introduction 
of dissipation or friction is expected to drive the system to 
simple fixed points in the phase space. This rather gen- 
eral expectation is surly true if the system phase space 
is bounded. However a much richer situation appears 
in the case of unbounded space, where unexpectedly a 
strange attractor can be induced by dissipation in an 
originally conservative system. To investigate this sit- 
uation we study the dynamics of particles on the Gal- 
ton board in the presence of constant external fields and 
friction. This board, introduced by Galton in 1889 ||, 
represents a triangular lattice of rigid discs with which 
particles collide elastically. For the case of free particle 
motion, the collisions with discs make the dynamics com- 
pletely chaotic on the energy surface as it was shown by 
Sinai [Q. In this paper we study how the dynamics of 
a charged particle in a presence of electric and magnetic 
fields is affected by a friction force Ff = —7V directed 
against particle velocity v (see Fig. 1). Without discs 
an external in-plane electric field E and a perpendicular 
magnetic field B create a stationary particle flow with 
the velocity Vf = (f + Fl)/7- Here f — eE is the effec- 
tive force, Fl = evf xB is the Lorentz force and e, m are 
the particle charge and mass. All perturbations decay to 
this flow with a rate proportional to 7 so that this lam- 
inar flow can be considered as a simple attractor. The 
effects of friction inside one cell of the Galton board at 
B = have been studied in || and it has been found that 
friction leads to the appearance of a nontrivial strange at- 
tractor. At present the effects of energy dissipation are 
actively investigated with the aim to construct equilib- 
rium and noncquilibrium steady states in a deterministic 
way (see J(|[jJ and Refs. therein). Here the Nose- Hoover 
thermostat with a momentum-dependent friction coeffi- 



cient allows to reach a number of interesting results with 
applications to molecular dynamics and nonequilibrium 
liquids ||. In our studies, contrary to 0], we concentrate 
mainly on the spatial structure of the turbulent chaotic 
particle flow appearing in the presence of friction. We 
show that the flow direction can be efficiently affected by 
a magnetic field. The obtained results describe the elec- 
tron dynamics in antidot superlattice which has been ex- 
perimentally realized in semiconductor heterostructures 
H . In such structures the effects of classical chaos play 
an important role ||] and the effects of friction we discuss 
here can appear for relatively strong electric fields. 




FIG. 1. Chaotic flow on the Galton board. Here the dis- 
tance between discs is R = 2.24, f = (-0.5, -0.5), B — 2 and 
7 = 0.1. Initially 200 particles are distributed homogeneously 
along a straight line segment in the upper right corner, their 
color homogeneously changes from red to green along this 
segment. 

To study numerically the dynamics of this model, we 
fix the disc radius a = 1 , and e = m — 1 so that the sys- 
tem is characterized only by the distance R between discs 
packed into equilateral triangles all over the (x, y) plane, 
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the friction coefficient 7, the external force of strength 
/ = |f| and the cyclotron frequency ui c — B jn|. For 
R < R c = the re is no straight line, by which a 

particle can cross the whole plane without collisions (at 
/ = 0, B = 0,7 = 0), and we start the discussion from 
this case. The particle dynamics is simulated numeri- 
cally by using the exact solution of Newton equations be- 
tween collisions, and by determining the collision points 
with the Newton algorithm. This way the trajectories 
are computed with high precision, and for example for 
7 = the total energy is conserved with a relative preci- 
sion of 10~ 14 |nj. A typical example of the chaotic flow 
formed by an ensemble of particle trajectories is shown 
in Fig. 1 . To illustrate the mixing properties of the flow 
we attributed a color to each trajectory that allows to 
follow their interpenctration and spreading. This figure 
shows that there is a certain penetration depth of one 
color into another, however this depth is finite since on 
average the initial color repartition is still visible. 
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FIG. 2. Example of a single trajectory for the case of 
Fig. 1 shown on small (main figure) and large (lower inset) 
scales. The drift velocity of the flow Vf « 0.13 is directed at 
angle a « 0.48 to f. Upper inset shows a single trajectory 
in the region (-160 < x < -140,-190 < y < -170) for 
7 = 0.004 with v f w 0.05. 

The properties of color penetration can be understood 
from the analysis of single trajectory dynamics. Such 
a typical example is presented in Fig. 2. It shows that 
the particle moves with an average constant velocity v/ 
under some angle a to the external force f, except for 
B = where a = 0. This drift velocity is constant only 
on average since on a smaller scale the particle moves 
chaotically between scatters following a strange chaotic 
attractor. In Fig. 2 for 7 = 0.1 the drift velocity is rela- 
tively large and the particle does not have enough time 
to move around many scatters in the direction perpen- 
dicular to the flow. As a result the penetration depth for 
color mixing is not very large. Surprisingly, for smaller 
friction, the drift velocity becomes smaller and the pene- 



tration depth increases so that the particle makes many 
turns around discs as it is shown in Fig. 2 for 7 = 0.004. 
This dependence is opposite to the case without scatters 
where vt drops with the 7 growth. 

This result can be understood on the basis of the 
following physical arguments (for B = see also ||). 
In the regime with weak friction the particles start to 
diffuse among discs in a chaotic manner with the dif- 
fusion rate D = vl/2 where v is the particle velocity 
and / is the mean free path |0|. For R ~ 1 we have 
I ~ R ~ 1 while the dependence of / on R will be dis- 
cussed in more detail latter. During the dissipative time 
scale r 7 = m/7 this diffusion leads to the particle dis- 
placement Ar ~ y/ Dm/ 7 along the direction of the drift 
velocity v f . This gives the change in the potential en- 
ergy U ~ /Arcoso; ~ f e f yj Dm / 7 where f e f = / cos a. 
In the stationary regime at time ( » t,, this potential 
energy should be comparable with the kinetic energy of 
the particle so that U ~ tow 2 , where v 2 —< v 2 + v 2 > is 
the average velocity square. Hence 



v 2 ~ (/cosa) 4/3 (Z/7m) 2 / 3 



(1) 



This relation allows to determine the drift velocity of 
the flow Vf. Indeed during the time r c between col- 
lisions the particle is accelerated by average forces f 
and Fl = eVf x B that gives the average drift velocity 
Vf = (f + Fl)t c /to. Since the dynamics is chaotic the 
direction of velocity is changed randomly after each col- 
lision so that Vf is accumulated only between collisions. 
The time r c is determined by the mean free path I and 
the average velocity v: t c = l/v = 2D/v 2 . Thus for the 
angle a between Vf and f we obtain the relation: 



tana = eBT c /m 



eB l 2 ^ 7 V3 



(2) 



(to/ cos a) 2 / 3 

The amplitude of fluctuations around this direction is 
Ar ~ yl Dm/ '7 which also determines the color mixing 
depth (see Fig. 1). 

From (pi) and (ph we obtain the drift velocity amplitude 



«/ = / e /r c /TO^? 2 / 3 (7/ e/ /m 2 ) 1 / 3 



(3) 



with f e f — / cos a. For B = the particles flow in the 
f direction. In this case their mobility is /1 = Vf/f — 
T c /m = D/{mv 2 /2). This is in fact the Einstein rela- 
tion according to which the mobility is given by ratio of 
the diffusion rate to the average kinetic energy (temper- 
ature) Jl|]. At B = the relations (§), @ and (§) are 
in the agreement with those in |^] and with the numer- 
ical data shown in Fig. 3. The values of Vf and v 2 are 
obtained from one very long trajectory (with a length of 
10 3 — 10 4 R ) or 10 shorter trajectories. Within statistical 
fluctuations this gives the same Vf and v 2 independent 
of their initial values. 

The relations ([!]), (||), ([|) allow to estimate the value 
of the Lyapunov exponent A. Indeed, the particle moves 
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with a typical velocity v and as in the case of the Sinai 
billiard with i ~ ~ 1 we have A ~ v/l. Therefore in 
the regime when: 



7 < 7c - yjmUj/l , (4) 

the value of A is much larger than the dissipation rate 
7/m. As a result for 7 <C 7 C the strange attractor is 
fat and its fractal dimension is close to the maximal di- 
mension 4, which is determined by the number of degrees 
of freedom (we remind that contrary to the nondissipa- 
tive case the energy is not conserved). For 7 ^> "f c the 
dissipation time r 7 becomes much shorter than the time 
between collisions r c . In this case the dissipation dom- 
inates chaos and the strange attractor degenerates into 
a simple attractor. For 7 < 7 C our numerical simula- 
tions performed with high computer accuracy show that 
trajectories remain chaotic for displacements from the 
origin being larger than 10 5 i?. Also at B = it can 
be shown that 7 C ~ mv/l ~ y/ fma/R for R 3> a and 
7c ~ y/mf/a 3 AR for AR = R - 2a < a. 



note that according to our data (see Fig. 4) the strange 
attractor exists even in the case R > R c — 4/V3 when at 
/ = 0, 7 = there are straight trajectories crossing the 
whole plane without collision. Apparently the contribu- 
tion of these orbits is not significant if 7 > and if f is 
not directed along these lines. 
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FIG. 3. Dependences of v 2 and Vf (inset) on / and 7 for 
R = 2.24, f/jf| = (-1,-1), B = and 0.001 < 7 < 0.4; 
0-5 < / < 32: circles show numerical data and lines show the 
slopes from (jl|) (main figure) and ^ (inset). 

The above changes in the mean free path I at B = 
also affect the drift velocity of the flow through the re- 
lation d3). Indeed, for AR < 1 we expect I ~ AR that 
gives Vf oc AR 1 . This dependence is close to the nu- 
merical data shown in Fig. 4 even if the numerical value of 
the exponent is approximated better by 0.5. In the other 
limit R 1, we have I ~ R 2 /a that gives Vf oc i? 4 / 3 . 
This power dependence is in a satisfactory agreement 
with the data in Fig. 4 although the numerical value of 
the exponent is closer to 1. We attribute these small de- 
viations in the exponent values to a restricted interval of 
variation in R. Actually we can not use very large/small 
values of AR since in these limits the value of 7 becomes 
comparable with j c and chaotic attractor disappears. We 



FIG. 4. Dependence of v f on AR = R-2 for 7 = 0.1, 5 = 
and f = (—0.5, —0.5): points give numerical data, lines show 
the slopes 0.5 and 1. 
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FIG. 5. Dependence of tana on scaled magnetic field 
B s = -B7 1/3 /r/ /3 horn (§) with f ef = /cos a for R = 4, 
f/|f| = (-1,-1) and 0.125 < B < 16, 0.14 < / < 5.65, 
0.002 < 7 < 0.1. Points give numerical data, line shows the 
average dependence tana = 1.6-B a - 

The introduction of magnetic field allows to change ef- 
ficiently the direction of the flow. The numerical data for 
the variation of tan a with the magnetic field and other 
system parameters are presented in Fig. 5. The aver- 
age dependence is in a good agreement with the equa- 
tion (||) for a large region of parameter variation where 
tana changes by two orders of magnitude. At the same 
time for moderate angles a < 1 the flow velocity Vf is 
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weakly affected by B. For example, for the case of Fig. 2, 
Vf remains practically the same for B = 2 and B = 
(vf w 0.15). For a > 1 magnetic field starts to change 
significantly v 2 and Vf in agreement with (Q) and (^). 
While the data in Fig. 5 in average follow the depen- 
dence (||) the fluctuations around the average are rather 
large. Their origin becomes clear from Fig. 6 where only 
the magnetic field is changed. In this case, tana has 
a pronounced peak which is located near the value of 
B where the cyclotron radius r c = v/B is equal to the 
disc radius. Indeed for r c > a a trajectory can make 
a full turn around a disc that allows to increase a and 
to reach a strong deviation of the global flow from the 
direction of the electric field. The growth of a leads to 
a drop of currant (conductivity) in the f direction and 
hence to the increase of resistivity. In fact, the peaks in 
resistivity near r c sa a had been observed experimentally 
H and explained theoretically |J in the linear response 
regime. Our data show that the peaks should also exist in 
the regime with strange chaotic attractor and relatively 
strong electric field where the I — V characteristics be- 
come nonlinear. For a denser package of discs the above 
peak is still present (R = 2.24 in Fig. 6) but it is less 
pronounced. It is interesting to note that in this case, 
a can even be negative so that the particles flow against 
the average Lorentz force. We note that the possibility of 
such a flow had been discussed for systems with Hamilto- 
nian dynamical chaos in |Q . Generally for a dense disc 
package the contribution of resonant orbits with r c « a 
starts to be significant and deviations from the average 
dependence (0) become rather strong. 
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FIG. 6. Dependence of tana on B for f = (—0.5,-0.5), 
7 = 0.03, R = 4 (points) and R = 2.24 (circles). Cyclotron 
and disc radiuses are equal at B ~ 1.6 where v 2 » 2.5. 

The above dynamical turbulent flow has rather inter- 
esting and unusual properties and it would be interesting 
to study it in experiments with antidot superlattices like 
in ||. The regime we discussed should appear when the 
steady state velocity v from (111) becomes larger than the 



Fermi velocity vp of the charge carriers. Thus, in addi- 
tion to (H), the condition v > vp determines the threshold 
for appearance of a strange attractor. The experimen- 
tal investigation of this phenomenon should also shade 
light on the role of quantum effects in such regime. We 
note that in the derivation of ([!]) - (||) we didn't use any 
specific properties of the disc distribution in the plane 
and therefore the results remain valid for randomly dis- 
tributed discs. Our study represents also certain interest 
for transport properties of neutral/charged particles sus- 
pended in a viscous flow streaming through a system of 
scatters. Indeed, a laminar stream with the velocity v s 
creates an effective force / = v s ^/ e f where j e f is the effec- 
tive friction created by the viscosity of the liquid. Such 
kind of transport can be studied experimentally with vis- 
cous liquids and its investigation can contribute to a bet- 
ter understanding of the interplay between dissipation, 
turbulence and chaos. 
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